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We consider a layered superconductor with a complex order parameter whose phase switches sign 
from one layer to the next. This system is shown to exhibit gapless superconductivity for sufficiently 
large interlayer pairing or interlayer hopping. In addition, this description is consistent with experi- 
ments finding signals of time reversal symmetry breaking in high-temperature superconductors only 
at the surface and not in the sample bulk. 



PACS numbers: 74.20.Kk, 74.50. -f-r, 74.72.-h 

It has been proposed that the many-particle ground 
state of the high-temperature superconductors breaks 
time-reversal symmetry T ||l|,^. As yet, no clear ex- 
perimental signatures of T-violation in the supercon- 
ducting bulk have been found H,^. Studies of T- 
violating ground states for model magnetic Hamiltoni- 
ans predict the nontrivially complex order parameter 
Ak = Ao(cos kx — cos ky + ie sin k^ sin ky) which is known 
as a d^2_y2+iedxy (d+id) order parameter This order 
parameter does not vanish anywhere on the fermi surface 
and therefore appears to be ruled out by experimental 
studies of high-temperature superconductors which in- 
dicate the presence of isolated nodes on the fermi sur- 
face [||-||. A different set of experiments finds signs of 
a nodeless T-violating order parameter at the surface of 
the sample. Evidence of a superconducting gap has been 
found at isolated point on the sample surface in STM 
studies H]. A gap feature of approximately 5 meV is 
always observed in c-axis tunneling between YBCO and 
a conventional superconductor Q. Anomalous distribu- 
tions of magnetic flux have been observed at crystal grain 
boundaries in YBa2Cu307 [|oj. This experiment is inter- 
preted as evidence for T-violation |jri|. In this Letter, we 
study the properties of a complex order parameter with 
a two layer unit cell. The phase of the order parameter 
changes sign from layer to layer. We show that inter- 
layer pairing and tunneling render this superconductor 
gapless at isolated nodes. Our calculation only describes 
bulk properties because it relies on the symmetries of the 
unit cell. Because these symmetries need not be satisfied 
at sample boundaries, our model is consistent not only 
with gapless superconductivity in the bulk but also with 
T-violating gapped superconductivity at the boundaries. 

Let us now investigate the behavior of multilayer 
T-violating superconductivity. Studies of coupled T- 
violating systems have found that it is energetically fa- 
vorable for the sign of the T-violation to be opposite 
in the two systems ||l2[. We therefore choose the order 
parameter to be Ak in the first layer and AJ^ in the sec- 
ond layer. We also assume an in-planc spectrum for the 
non-interacting electrons of ek, a hopping matrix element 
between layers of tj^k, and a bilayer order parameter of 



A^k- The standard approach is to define quasiparticle 
creation and annihilation operators a^„^ (i = 1,2 for the 
first and second planes respectively, a =T:i) [|l3|. The 
Hamiltonian is then 

kcri 
kcr 

- X [A_Lk(ak|ia^^j^2 + a-ki2akTi) + C-C 

k 

Ak(a[|iaL^j^;^ -I- a_ki2akT2) + c.c . (1) 

This expression can be rewritten using Nambu notation 
as 
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(2) 



where 



Qk = 



r 


-Ak 


t±k — Aj^k " 


-a;; 


-a 


-A^k -t±k 


^Ik 


-A_Lk 


Ck -A* 


L-Alk 


-^Ik 


-Ak -Ck . 



OkTi 
a-kii 

OkT2 

t 

ki2 



a 



(3) 

$k = Ck — M = — 2i(cos kx + cos ky) — fi, Ak = Ao(cos k^ — 
cos ky + ie sin k^ sin ky ) , and both A j^k and t j^k are taken 
to be real. Writing Ak = A^k + «A/k, we find that Qk 
has energy eigenvalues E± (k) that satisfy 



Aik + 4k 



± 2^(aUk + AikAijk)2 + A2^(i^i^ + A^, 
The spectrum for i?_ (k) vanishes when 

Ck — i±kak , Aflk = A_Lkak, 

and 
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which is consistent with experiments that find a bulk or- 
der parameter with a d-wave like node P-|^ JT^ , |l5t . If the 
phase of the order parameter were the same in neighbor- 
ing planes, the corresponding energy eigenvalue i?±(k) 
would satisfy 

E'lik) = (a ± t±k)' + (A/k)' + (Aflk T A^k)" (7) 

instead of Eq. (^) and hence be nodeless except for a 
special value of A^k- 

In order to have a nodeless order parameter, we must 
be able to satisfy Eq. @-(^). The constraint that be 
real requires 



^lk + 



^J^k ^ ^7k 



(8) 



This intuitive result tells us that the combined effect of 
interlayer pairing and interlayer tunneling has to be suf- 
ficiently large in order to overcome the gap due to the 
imaginary component of the order parameter. The exis- 
tence of a solution to Eq. (||) also requires Ak > Aj^k- 

The zeros of Eq. (^) need not lie on the 45° {k^ — ±ky) 
nodal line of the pure (1^2 _y2 order parameter. To first 
order in t^/t and A^/Aq and in the t » eAo limit, the 
nodal wavevector k lies at 
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The node is shifted from the 45° nodal line to an angle 

akAj_k 



a,TCtSin{ky / kx) w 7r/4 + 
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(11) 



which is independent of t_Lk to lowest order. In contrast, 
the magnitude of k. 



|k|«V2(| + ^^ 



(12) 



is independent of Aj^k to lowest order. There are now 
two nodes because may be either positive or negative. 
These effects (node splitting and angular shift of 10°) 
have been observed in angular resolved photoemission 
studies of Bi2Sr2CaCu208 ||l6[. A similar analysis of the 
momentum shifts has been performed for the case of a 
purely real order parameter by P. A. Lee and coworkers 
||17| , |l8| . Our expressions are valid in this case with ak = 
±1 in Eq. (||)-(|l^). Therefore, the bulk photoemission 
experiments cannot distinguish qualitatively between a 
purely real and an alternating complex order parameter. 

We may also consider other alternating complex order 
parameters. For example, we might consider an alternat- 
ing order parameter that is purely complex. Requiring 



that Ai{k = is equivalent to imposing an additional 
constraint upon Eq. (|^). This equation will only have 
a solution in the unlikely case that either Aj^k = or 
Q!k — 0. An order parameter with s + id symmetry 
^ would only satisfy Eq. (|)-(|) if the s component 
is smaller than Aj^k- 

The parameters we choose depend on the material we 
wish to study. For example, neighboring planes in YBCO 
may be so strongly coupled that they prefer their order 
parameters to have the same complex phase. In that case, 
the bi-layer in YBCO may be treated as a single compos- 
ite T-violating superconducting layer coupled to another 
composite layer with the opposite sign of T-violation. In 
order to observe a node, Eq. (H) must be satisfied for 
the coupling of the composite layers. This description 
should be consistent with the spin gap calculations for 
bilayer systems One may also consider 3-layer ma- 
terials (which arise in the Bi, Tl, and Hg cuprates) or 
even the effect of intermediate normal layers sandwiched 
between superconducting layers with alternate signs of 
T-violation. 

It is expected that a T-violating layer will couple to 
an electron or neutron spin as if it induced a local mag- 
netic field. At long wavelengths, this effect would be 
suppressed because of the alternating chirality. However, 
a neutron with momentum q = (0,0, tt/c?), where d is 
the distance between planes, will couple to the chirality 
oscillation. If the coupling is strong enough, a peak at 
low energies near this momentum should be observable in 
neutron scattering. The size of this feature will depend 
on the coupling strength, which has not been accurately 
calculated. However, muon spin rotation experiments, 
which probe the local magnetic field, fail to find an effec- 
tive magnetic field larger than that due to the Cu nuclei 
[Q. In contrast, similar studies of UPt3, which is believed 
to have a d + id order parameter that does not alternate, 
saw this effect 

Our postulated state can have interesting effects in ex- 
ternal magnetic fields if the energy scales allow these ef- 
fects to occur below Hc2- The energy gained aligning 
the chirality in all of the layers with a c-axis magnetic 
field may be contrasted with the energy gained alternat- 
ing the chirality in neighboring layers |12| . If the cou- 
pling to the magnetic field is unusually large, there may 
be an H* < Hc2 above which it becomes energetically 
favorable for the layers to violate T without alternat- 
ing phase. Hysteretic effects would also be expected for 
samples cooled in sufficiently large magnetic fields. 

There are three modes made from linear combinations 
of the phases of the order parameters Ak (the first layer), 
A^ (the second layer), and Aj^k- The mode where the 
phases are added is the charged Anderson-Bogoliubov 
mode which is lifted to the plasmon frequency. The other 
two linear combinations will have interesting dynamics 
and, as recently pointed out by Kuboki et al. jl^, should 
be observable as a Raman signal. Investigation of these 
modes may provide a criterion for differentiating between 
our oscillating d + id order parameter and the purely real 
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The effect of a 



one proposed by K. Kuboki et ah 
complex order parameter on Josephson tunneling 24-26|] 
should also be studied. 

The T-violating component of the alternating order 
parameter should have a transition temperature T* < 
where is the bulk transition temperature. The ex- 
perimental results that we attribute to T-violation are 
therefore expected to disappear for temperatures above 
T* . In fact, recent experiments in oxygenated samples 
of YBCO have foimd a second transition at 30K in T2 
measurements p7| ]. 

In this Letter we propose that a certain class of T- 
violating superconducting ground states is consistent 
with numerous experiments. There remain significant 
inconsistencies. Most notably, several experiments have 
been interpreted as finding zero Hall conductivity. In 
contrast, T-violation is consistent with a nonzero value 
of the Hall conductivity Uxy It is difficult to compare the 
limits on the Hall conductivity set by experiment with 
theoretical predictions because this quantity has not yet 
been accurately calculated ^. Gauge field fluctuations 
are expected to reduce the size of the hall conductivity 
P, p^|j2q | . It is also known that fluctuation effects tend to 
reduce the quasiparticle gap in single layer anyon models 
pof and hence the size of the T-violating order parameter 
A/. This effect is expected to be enhanced in a double 
layer anyon model p3[ . 
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